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Abstract 

The adiabatic sound velocity and compressibility for harmonically trapped 
ideal anyons in arbitrary dimensions are calculated within Haldane frac- 
tional exclusion statistics. The corresponding low-temperature and high- 
temperature behaviors are studied in detail. To compare with the experi- 
mental result of unitary fermions, the sound velocity for anyons in the cigar- 
shaped trap is derived. The sound velocity for anyons in the disk-shaped 
trap is also calculated. With the parameter g = 0.287, the sound velocity 
of unitary fermions in the cigar-shaped trap modeled by anyons is in good 
agreement with the experimental result, while that of unitary fermions in the 
disk-shaped trap is Vq/vf = 0.406 with Fermi velocity vp- 
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1. Introduction 

The statistical behaviors of a quantum many-body system can be suc- 
cessfully described by Bose-Einstein or Fermi-Dirac statistics. Beside these 
two statistics, a new quantum statistics called Haldane fractional exclusion 
statistics had also been proposed by Haldane 

Following the distribution function of the Haldane anyon gas given by Wu 
jif , the thermodynamics of ideal anyons was a focus of theoretical attention 
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in the past few years 0, S 0, S, H B i, 0, Q H [3 Q E, M, [l7 
For example, the thermodynamics of the ideal anyons in two dimensions has 
been investigated [s], 0, [sj, and the analytical expressions for the chemical 
potential, internal energy, entropy, isochore and isobar heat capacities, and 
Joule-Thomson coefficient of a three-dimensional ideal anyon gas have been 
derived in the previous works 0, 0, 0, [oj. Furthermore, the thermodynamics 
for a homogeneous ideal gas within fractional exclusion statistics has been 
studied in some detail in arbitrary dimensions by Refs. 10|, lUl, ll2|, ll3|, uM, 



15l . Il6| . For a harmonically trapped d-dimensional ideal anyon system, the 



internal energ y a nd isochore heat capacity as a function of temperature were 



given in Ref. [17 



In recent years, the strongly interacting physics in ultracold fermions has 
attracted much attention both experimentally and theoretically 18|, ll9|, |20|, 



2ll . I22I I23I . l24j . I25I . l26l . 1271 ] . The ability to widely tune the effective interaction 
between two fermions by a broad Feshbach resonance in the atoms of ^Li 
and ^"K has permitted the experimental observation of a smooth evolution 
of the Fermi gas from the Bardeen-Cooper-Schrieffer (BCS) weakly attractive 
regime to a molecular Bose-Einstein condensate (BEC) regime. On the cusp 
of the BCS-BEC crossover, there is a strongly interacting regime that is the 



Fermi gas [20 



so-called unitarity limit regime [18|, |19|, which corresponds to the unitary 

the 
The 
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According to the theoretical universal hypothesis 
unitary Fermi gas can exhibit a universal thermodynamic behavior, 
universality shows that the ground-state energy /ig of a homogeneous gas at 
zero temperature should be proportional to the free Fermi energy ep with 
a proportionality constant ^ (C=l + P) which is the universal many-body 
constant. The ratio of the ground-state energy to the Fermi energy for the 
unitary system trapped in a harmonic trap is ^^^^ with the density functional 

$,^^'^eF, where ep is the trapped ideal Fermi energy. 



approach [22 



Thus, /io 

the thermodynamics of a trapped unitary Fermi gas has been 



Meanwhile 

experimentally measured 
velocity of a unitary Fermi gas inside a strongly elongated cigar-shaped trap 



23l . I24J . I25I |26| . Particularly, the ultracold sound 



had been measured by the Duke group [25|, |26|. These recent experiments 
on the thermodynamics of ultracold Fermi gases provide an evidence for 
the universality in the laboratory. Furthermore, the important universal 
parameter 1^ (or P) can be extracted from these experimental techniques and 
measured physical quantities. In the measurement of sound velocity in a 
unitary Fermi gas, Joseph et al. gave (3 = —0.565 ± 0.015 [25 . 



Physically, the thermodynamic properties of a unitary Fermi gas are be- 
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tween those of Fermi and Bose gases 27|. Meanwhile, the thermodynamic 
behaviors of the anyon gas obeying the Haldane fractional exclusion statis- 
tics are also between the fermionic and bosonic ones Due to the scale 
invariance, the thermodynamic quantities of the unitary Fermi gas and the 
Haldane anyon gas are both related with the corresponding quantities of the 
ideal Fermi gas at zero temperature. Therefore, it is assumed that the three- 
dimensional anyon gas obeying fractional exclusion statistics can be used to 
model the statistical behavior of a Fermi system at unitarity 
In addition, the finite-temperature internal energy and entropy of the ideal 
anyons within Haldane statistics are in good agreement with experimental 
data of the unitary fermions for a given statistical parameter g f^. Here, 
the value of g is determined by the ground-state energy and the value of the 
universal many-body constant ^ extracted from the ultracold experiment. 
In thermodynamics, the adiabatic sound velocity and compressibility for 



the homogeneous quantum gases are two important quantities [ll|, |28|, l29|, 



30, |3l|, |32|. The corresponding results of the ideal Fermi and Bose gases were 



obtained by Ref . |28[ . The adiabatic sound velocity and compressibility of a 
statistically interacting quantum gas in a generalized nonlinear Schrodinger 
model were calculated in Ref. [29|]. More recently, for a strongly interacting 
Fermi gas, the adiabatic sound velocity was discussed within the quasilinear 



approximation framework [30|, |3l| , and it was also calculated with an effective 



field theory near the Feshbach resonance regime [32 



In this Letter, the analytical expressions of the adiabatic sound velocity 
and compressibility for the harmonically trapped ideal anyons in arbitrary 
dimensions are derived within the Haldane fractional exclusion statistics. 
The corresponding low-temperature and high-temperature expansions of the 
thermodynamics are given. We also investigate the zero-temperature sound 
velocities for the three-dimensional anyon gases in the cigar-shaped and disk- 
shaped harmonic traps to model the corresponding physical quantities for a 
real unitary Fermi gas. 

The outline is as follows. In Section |2l the grand canonical partition 
function of the Haldane fractional exclusion statistics is given. The adia- 
batic sound velocity and compressibility of a li-dimensional trapped ideal 
anyon gas are derived analytically in Section [31 The corresponding numer- 
ical results are also evaluated in this section. To compare with the recent 
experimental result of ultracold fermions at unitarity, the zero-temperature 
sound velocity for a three-dimensional anyon system in a strongly elongated 
cigar-shaped harmonic trap is derived in Section |H In Section the sound 
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velocity for a three-dimensional anyon gas in a disk-shaped trap is also cal- 
culated analytically. Summary and conclusions are presented in Section El 



2. Grand canonical partition function 

In the fractional exclusion statistics, the number of microscopic quantum 
states of identical particles occupying a group of G states is |2[ 

iim[G,-gN,-{l-g)]V ^ ^ 

where the parameter g connects the change in the number of available states 
when one particle is added. 

In order to find the grand canonical partition function, let us first focus 
on the average occupation number (ra) = Ni/Gi. According to the Lagrange 
multiplier method 5\aW — a6N — P6E = 0, the average occupation number 
(n) is derived as Q 

where u and g satisfy the relation e = ksT In [uj^{1 + cj)^~^] + with the 
single-particle energy e, chemical potential fi, temperature T, and Boltz- 
mann constant ks- a = —fi/{kBT) and (3 = 1/(/cbT) are the two Lagrange 
multipliers. = ^ ■ Ni is the total particle number and E = ^ . e^A^j is the 
total energy. One has g = for bosons and g = 1 for fermions. Analogous to 
the Fermi-Dirac statistics, the zero-temperature average occupation number 
is (n) = with e > /i, and (n) = l/g with e < /i. 

Setting the grand canonical partition function as S = TT^S p-, the mean 
occupation number (n) of momentum p turns out to be [l6l . |28 

, ^ 1 fd\nE\ 

(n) 



(u) + g) duj 
where the fugacity z is defined as 



V / z,T,all other p 

uj{l + uj)dlnEp 



= exp(^) =u:,^{l + u;oy-\ (4) 
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Combining Eqs. ([2]) and (E]), one can conclude that 

2p = l + -, (5) 

where the integral constant is set to be vanishing. Correspondingly, the grand 
canonical partition function S can be identified to be 



3. Adiabatic sound velocity and compressibility 

We now discuss the trapped system thermodynamics. The geometric 
mean of the trap frequencies is defined as w = iJXi^i)^^'^ = 1;2, ■ ■ -^d). 
The density of states is 

where h = h/{2'K) is the reduced Planck constant, and rfoQis the gamma 



function. The corresponding system volume isV = w~ |17l.l33|. Notice that 
the volume here is no longer the usual thermodynamic variable for a system 
in rigid walls without the external potential. It is called the harmonic volume 



33|. We will discuss the physical meaning of the harmonic volume in the 
following of Eq. (fT3|). 

Utilizing Eq. (jH]) and integrating by parts, the logarithm of the grand 
canonical partition function S can be represented by turning the sum of 
quantum state into integral 



oo 



InS =2 / D(e)ln( l + -\de 







d 



where 



2(^) Ga-,M9). (8) 



1 r ^'~'dx 

Gd{,z,g) = —— ■ (9) 

V[d) Jo ^ + 9 
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is the Calogero-Sutherland integral function 
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It satisfies the relation 



dGd{z,g) _ \ Gd-i{z,g), if ^ 2, 
^ dz 1 if d = 1. ^ ' 

Further, the grand thermodynamic potential is = —PV = — /^^TlnS, 
where P is the press. The expressions for the particle number N, and internal 
energy E are 



N =z 



d\nE 

dz 



T,V 

d 



E = 



2(^) Gd{z,g), (11) 
91nH 



dp 



d 



= 2t/(^-^J kBTGd+i{z,g). (12) 

From the thermodynamic formula E — TS + PV = N/j,, one can obtain the 
entropy per particle 



— = Kr 



{d+l)Gd+i{z,g) 
Gdiz,g) 



(13) 



Before making further analysis, let us discuss the physical meaning of the 
harmonic volume V = w^'^. Since the temperature T and chemical potential 
/i are intensive, it is found that ro"'' must be extensive from the formulas of 
A^, E and S. At a given temperature, f2, iV, E and S are all proportional 
to w~'^. Thus, the harmonic volume V = w~'^ does make physical sense as a 
system volume" since small frequencies of trap imply large actual volumes 



33| . As shown in Eq. fllip . Nw'^ — t- constant in the thermodynamic limit, 
A^ — 7- oo and — t- 0. Therefore, the particle number density n is set to be 
n = Nvj'^ which will be used in the derivation of the adiabatic sound velocity 
below. 

In terms of hydrodynamics, the expression for the sound velocity is 
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with the particle mass m. 

According to the Jacobian techniques in thermodynamics, one has 



djPV) \ f d{S/N) \ 

\dn)~\ dN J^y f d{s/N) \ 



(15) 

N,V 



With Eqs. (dH), ([IS]), dH]), and ([I5D, we finally obtain the 

normalized sound velocity 



VF V 2d Tf GAz,g) ' ^ ' 



where vf = \/2ep/m is the zero-temperature Fermi sound velocity of a 
harmonically trapped ideal Fermi gas. 

As indicated in Eq. f[T^ . we have transformed the intensive thermo- 
dynamic quantizes P and n into the extensive quantizes VL = —PV and 
= nw~'^, which are not localized. That is to say, Q and A^ are indepen- 
dent of location. Essentially, the sound velocity obtained in Eq. ( IT6!) is an 
average velocity across the harmonic trap. Thus, v is independent of the 
position in the trap. 

Similar to the derivation of the sound velocity, the adiabatic compress- 
ibility turns out to be 



- _1 
~ ~V \dP 



s 

d 1 Gd{z,g) 
d+ InksT Gd+i{z,g)' 



;i7) 



3.1. Low-temperature behaviors 

In order to derive the low-temperature expressions of the thermodynamic 
quantities, we consider that the zero-temperature particle number is equal 
to the finite-temperature one. 

At zero temperature, it is worthy noting that the particle number can be 
alternatively given by 

N=- D(e)de = , 'f, -, (18) 
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where ?p obeys the relation ?p = g^^^ep with the ideal Fermi energy ep 
[T{d + l)A^/2]^/'^ficj7 in a harmonic oscillator. 

By substituting Eq. f|T8|) into Eq. f|TT]) . one gets 



T 
7> 



r(c/ + i)^ 



(19) 



with the Fermi characteristic temperature Tp = [T{d + l)N/2Y^'^hzu/kB for 
a trapped ideal Fermi gas. 
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At very low temperature, the g) can be expanded as [ll|, ll2|, |13|, Il4j . 



(In z)" 

^?r(n + i) 



1 + 



TT^ (7n(n — 1) 
~6 {Inzy 



+ 



(20) 



It should be pointed out that Gi{z,g) can not be expanded by applying 
Eq. fl2U]) . In fact, its expression with uq or temperature T can be directly 
given by Eqs. ([9]) and 



GAz, 



In 



1 + CJQ 



Tp 
T ' 



(21) 



Therefore, substituting Eq. ( l20i) into ( |T9l) and ( |T2l) . one obtains the low- 
temperature analytical expressions of the chemical potential and internal 
energy as follows: 



ep 



i^tl^W l-2/d 

6 ^ 



T 
7> 



+ ■ 



(22) 



E 



d 



Nep d + l 



l/d 



1 , (^-lK% i-2/d 



T 
7> 



+ 



(23) 



For the special case c? = 1, the dominant asymptotic expansions in Eqs. 
and fl23|) vanish. Therefore, we obtain the very weak results ^/ep ~ g and 
E/{Nep) ~ ^/2, as T ^ 0. 

By applying Eqs. (EO]) and (El]) to Eqs. (JTH]), ([161), and dlT]), one has 



S 



l-l/d 



7> 



+ 



(24) 



V 

Vf 



'2d 



T 



(25) 



d 



ng^/'^eF 



6 



+ ■ 



(26) 



3.2. High-temperature behaviors 

At very high temperature, the Calogero-Sutherland integral functions 



have the power series expansions as follows UJ, Il2 

-v2 



Therefore, one gets 



z + -[l-2g) + 



(27) 



V 

Vf 



d+l T 
2d 7> 



1 + ^(2^?-!) + 



(28) 



d 



d + l nksT 



(29) 



From Eqs. fl28l) and fl29l) . it is found that f/f^ — ^ cxo and Kg — )■ when 
T/Tp -> 00. 

5.5. Numerical results and discussion 

Now, we will present the numerical results in this subsection. 

3.3.1. Adiabatic sound velocity 

The normalized adiabatic sound velocity for a trapped ideal anyon system 
versus the rescaled temperature can be calculated from Eqs. (IT^ . ( IT^ . and 



As indicated by Fig. [H the adiabatic sound velocity of a trapped ideal 
anyon gas increases with the increase of the rescaled temperature. The 
curves for trapped anyons get closer to that of the ideal fermions in the 
high-temperature Boltzmann regime. The adiabatic sound velocity given by 
this model is lower than the ideal fermionic one, and it goes to a smaller 
value for a smaller value of g at the same temperature. 
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Figure 1: The normalized adiabatic sound velocity is plotted as a function of the 
rescaled temperature ((a) d = 3; (b) d = 2; (c) d = 1). The solid curve denotes that 
for the trapped ideal fermions, and the dashed curves denote the ones for trapped 
ideal anyons with statistical parameter g = 0.1,0.3,0.5,0.7,0.9, respectively. 



3.3.2. Adiabatic compressibility 

The adiabatic compressibility for a trapped ideal anyon system versus the 
rescaled temperature can be calculated from Eqs. ( IT71) . ( fT9l) . and ( l26l) . 

Fig. |2] shows that the adiabatic compressibility of a trapped ideal anyon 
gas decreases with the increase of temperature. In the high-temperature 
limit, the curves for trapped anyons get closer to that of the ideal fermions, 
and they all tend to zero. The adiabatic compressibility increases with the 
decrease of g for fixed T. 
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Figure 2: The normalized adiabatic compressibility versus the rescaled tempera- 
ture ((a) d = 3; (b) d = 2; (c) d = 1). The line styles are similar to the ones in 
Fig. [TJ It is indicated from Eq. (j26p that the zero-temperature compressibility is 
Kg = d/^nep) for a trapped ideal Fermi gas. 



4. Sound velocity for a three-dimensional anyon gas in a strongly 
elongated cigar-shaped trap 

Experimentally, the measurement of ultracold sound velocity for a unitary 
Fermi gas in a strongly elongated cigar-shaped trap had been given in Refs. 



25|,|26j. In order to compare with the experimental datum, let us consider the 
case of a three-dimensional anyon system in a strongly elongated cigar-shaped 
harmonic trap V{z,r±) = (m/2)(a;^z^+a;^r^) with u± ^ Uz- Here, r± is the 
radial axis, u± is the radial trap frequency, and Uz is the axial trap frequency 
along z axis. For making an analytical approach available, we assume the 
trapping potential is z independent, i.e., V{z,r±) = V{r±) = {m/2)uj\r']_. 
Within the local density approximation, the chemical potential of a cigar- 
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shaped harmonically trapped system in equilibrium can be expressed as 



and the radial Thomas-Fermi radius R± is determined by 

1 

2' 



±5 



(30) 



(31) 



where /iq is the zero-temperature chemical potential, and ^o{r±) is an effective 
chemical potential which depends on r^. 
With Eqs. ([20]) and (EH), one can get 



no(r^) 



nh 



1^0 



3/2 



1 - 



Ri 



3/2 



(32) 



where rih = (2m/io)^^^/(37r^^^) that includes two degrees of the spin degen- 
eracy, and R± = [2iXQ/{muj'j_)Y^'^. With the Thomas-Fermi approximation, 
the ground-state particle number Nq of the system is given by 



Nq = — 2-Kno{r±)r±dr_i_ / dz 
9 Jo J-Z/2 

_ Z8(2m)i/2 



5/2 

2 ■ 



(33) 



If the cigar-shaped trap is equivalent to a cylinder, the constant Z in Eq. 
fl55]) is the length of the cylinder. 

Therefore, one can derive the corresponding normalized sound velocity at 
zero temperature as 



Vf 



No duo 



vp\m ONq 



(34) 



For a harmonically trapped three-dimensional anyon gas at zero temperature, 
it is indicated from Eq. fl2^ that /ig = g^^^ep- Correspondingly, with Eq. 
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one gets 



As discussed in the third paragraph of the introduction part, there is /xq = 
^^^"^ep for a zero-temperature trapped unitary Fermi gas. Comparing /iq = 
g^^^ep with yUo = ^^^"^ep, one can find that the relation between the statistical 
parameter g and the universal constant ^ is g = ^^/^ (sl, 0, Ist- Consequently, 
the final result is 

— = ^> 36 



which is consistent with Eq. (2) in Ref. |25| based on the hydrodynamic 



approach that given by Capuzzi et al. |34 . 

The recent experimental measurement on ultracold sound velocity for a 
unitary Fermi gas in a strongly elongated cigar-shaped trap provided ^ = 
1 + (3 = (1 — 0.565) =F 0.015 (255. Therefore, one can determine the statistical 
parameter g = ^3/2 = 0.287. 



5. Sound velocity for a three-dimensional anyon gas in a disk- 
shaped trap 

For the anyon gas in a disk-shaped harmonic trap with Uz ^ uj±, we 
assume that V{z,r±) = V{z) = {m/2)uj1z^ . The chemical potential is /xq = 
jlo{z) + (m/2)(X'^-2^, and the axial Thomas- Fermi radius is defined as Rz = 
{l/ujz)\/2^,Q/m. 

The corresponding ground-state particle number takes the following form 



1 



R 



Nq = — I no(z)dz I 27rr±dr± 

ttR^ m 2 
" g 27r;i3^/°' 



(37) 



where 



no{z) = Uh 



z 

'r. 



3/2 



(38) 
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If the disk-shaped trap is regarded as a cyhnder, the constant 'nE? in Eq. 
( P7|) is the cross-sectional area of the cyhnder. 

Similarly, the reduced zero-temperature sound velocity is given by 

!^ = IM (39) 
Substituting yUo = g^^^^F and g = ^^^"^ into Eq. fl39|) . one obtains 

i^^ii^. (40) 

Vf 2 ^ ^ 

With g = 0.287, the normalized sound velocity of a unitary Fermi gas in a 
disk-shaped trap is vq/vf = 0.406. 



6. Summary and conclusions 

The analytical expressions of the adiabatic sound velocity and compress- 
ibility for a harmonically trapped ideal anyon gas in arbitrary dimensions are 
derived within the Haldane fractional exclusion statistics. The corresponding 
low-temperature and high-temperature expansions are also studied. 

With careful numerical study, we find that the adiabatic sound velocity of 
a trapped ideal anyon gas is a monotone increasing function of temperature 
for a given g, and the curves for the trapped anyon gases get closer to that of 
the ideal fermions in the high-temperature Boltzmann regime. The adiabatic 
sound velocity given by this model is lower than the ideal fermionic one, and 
it goes to a larger value for a larger value of g at the same temperature. 
However, the adiabatic compressibility is a monotone decreasing function of 
temperature for a fixed g, and it goes to a larger value for a smaller value of 
g at the same temperature. The adiabatic compressibility curves all tend to 
zero in the high temperature limit. 

The ultracold sound velocity for a three-dimensional unitary Fermi gas in 
a strongly elongated cigar-shaped trap modeled by anyon statistics is consis- 
tent with the experimental measurement for a given g. The sound velocity 
for a three-dimensional anyon gas in a disk-shaped trap is also obtained. 
With g = 0.287, the normalized sound velocity of a unitary Fermi gas in the 
disk-shaped trap is Vq/vf = 0.406. 
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